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ABSTRACT 

We explore what dominant physical mechanism sets the kinetic energy contained in neutral, atomic 
(H I) gas. Supernova (SN) explosions and magneto-rotational instability (MRI) have both been pro- 
posed to drive turbulence in gas disks and we compare the H I line widths predicted from turbulence 
driven by these mechanisms to direct observations in 11 disk galaxies. We use high-quality maps of 
the HI mass surface density and line width, obtained by the THINGS survey. We show that all sample 
galaxies exhibit a systematic radial decline in the H I line width, which appears to be a generic property 
of HI disks and also implies a radial decline in kinetic energy density of HI. At a galactocentric radius 
of - often comparable to the extent of significant star-formation - there is a characteristic value of 
the HI velocity dispersion of 10 ±2 km s _1 . Inside this radius, galaxies show HI line widths well above 
the thermal value (corresponding to ~ 8 km s~ x ) expected from a warm HI component, implying that 
turbulence drivers must be responsible for maintaining this line width. Therefore, we compare maps 
of HI kinetic energy to maps of the star formation rate (SFR) - a proxy for the SN rate - and to 
predictions for energy generated by MRI. We find a positive correlation between kinetic energy of HI 
and SFR; this correlation also holds at fixed Em, as expected if SNe were driving turbulence. For 
a given turbulence dissipation timescale we can estimate the energy input required to maintain the 
observed kinetic energy. The SN rate implied by the observed recent SFR is sufficient to maintain 
the observed velocity dispersion, if the SN feedback efficiency is at least £sn — 0.1 x (10 7 yr/r^), 
assuming td — 10 7 yr for the turbulence dissipation timescale. Beyond r25, this efficiency would have 
to increase to unrealistic values, e > 1, suggesting that mechanical energy input from young stellar 
populations does not supply most kinetic energy in outer disks. On the other hand, both thermal 
broadening and turbulence driven by MRI can plausibly produce the velocity dispersions and kinetic 
energies that we observe in this regime (> r2s)- 

Subject headings: galaxies: evolution - galaxies: ISM - galaxies: kinematics and dynamics - galaxies: 
dwarf - galaxies: spiral - stars: formation 



1. INTRODUCTION 

The extended HI disks of spiral galaxies exhibit 
characteristic 21 cm line widths of 5-15 km s _1 
on small scales. Line widths are observed to be 
~ 6-10 km s _1 even beyond r25, where the star 
formation rate (SFR) and the associated feedback 
are presumably inefficient (|van der Kruit fc Shostakl 
" 19821 iDickev. Hanson, fc Helqul 119901: iKamphuisI 119931 : 
Rownd. Dickev. fe Heloulll994 Ivan Zee fc Brvamm999f h 
In some observed cases, the HI disks exhibit an overall 
outward decrease of the velocity dispersion to 6-8 km s ~ 1 
(jBoulanger fc Viallefondlll992t iPetric fc Rupenll2007ft . 

The observed line widths are attributable in part to 
thermal broadening and in part to turbulence. At pres- 
sures typical of disk galaxies, interstellar atomic gas 
has two temperatures at which it can rem ain in sta- 
ble thermal equilibrium (IWolfire et al.l IT995) . The gas 



is observed to be distributed with roughly 60% of the 
mass evenl y distributed between the equilibrium tem- 
peratures dKulkarni fc Heilesl Il987t IDickev fc Lockman] 
1990; Dick ev fc Brinks! 1993h and 40% transiting between 
them (|Heilesll200lD . The cold (~ 100 K) neutral gas emits 
lines with a characteristic line width of ~ 1 km s^ 1 , while 
the warm (~ 8000 K) neutral gas has a line width of 
~ 8 km s _1 . Spectral lines observed to be broader than 
~ 8 km s^ 1 can be attributed to turbulence - stochas- 
tic fluctuations of velocity and pressure on scales much 
larger than the mean free path of the gas. Heating by a 
galactic interior or e xtragalactic b ackground of ultravi- 
olet (UV) radiation (jSch aye 2004|) may be significant in 
the outermost regions of H I disks, where velocity disper- 
sions approach thermal values for the warm phase. 

There are several energetic sources that can induce 
turbulent motions in the interstellar medium (ISM), 
among which feedback from recently formed stellar 
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iBalbus fc Hawlevl 119911 . 119981 : ISellwood fc Balbusl 119991 ) 
are some of the most important. The stirred tur- 
bulence drives motions at smaller and smaller scales 
until viscosity dissipates the motion into thermal en- 
ergy. This hierar chical description was formulated by 

iKolmogorovl (jl941) for an incompressible turbulent fluid, 
though it applies well to the small scale dynamics of the 
gas i n galaxy disks (jElmegreen. Elmegreen, fc Leitrierl 
2003). Ultimately, the ISM radiates away the ther- 
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mal energy, e.g, by line cooling. S ince the typical ra- 
diative cooling time (~ 10 3 years; Wolf ire et al] l2003h 
is shorter than the time scale for turbulent energy to 
decay (~ 10 Myr; IMac Low! I1999D , the turbulent en- 
ergy must be continuously replenished on the decay time 
scale. The direct identification of the energy sources 
that sustain the turbulence at this rate has drawn con- 
siderable interest for well ov er three decades, but has 
not yet been concluded (e.g. IMac Low fc Kless"enl 12004 
lElmegreen fc Scaloll2004 . 

Numerical simulations have demonstrated that both 
hydrodynamical and magnetohydrodynamic turbulence 
decay over a relatively rapid time scale, compara- 
ble to the loca l dynam ic al time, i.e. ~ 10 Myr 
(IMac Low et al] 119981; IStone. Ostriker. fc Gammid 



1998 1 ; IMac Low! Il999t iPadoan fc Nordlun 



1999: lAvila-Reese fc Vazquez-Semadenil 120011; 
Ostriker. Stone, fc Gammidl200l 7 thus requiring a driv- 



ing mechanism that continuously pr ovides energy to the 
medium to maintain a steady state (Gammie fc Ostriker! 
19961: IMac Lowlll999ft . 

Which mechanisms drive ISM turbul ence remains a 
matter of debate. As first envisioned bv ISpitzerl (|1978h . 
star formation itself injects energy back to the ISM 
and drives turbulence through stellar winds of O/B and 
Wolf-Rayet stars, UV ionizing radiation, and super- 
nova (SN) explosions (e.g. IChiang fc Prendergastj 
1985t iRosen. Bregman. fc Normanl 119931; 



Kim. Balsara. fc Mac Low! 120011; iDib. Bell, fc Burkertl 



Joung fc Mac Low! 120061 : lAvillez fc Breitschwerdtl 
The energy inpu t from SNe likely dom- 
inates stellar feedback (|Kornreich fc ScaioT 120001 ; 
IMac Low fc Klessenll2004T ). The overall galaxy kinemat- 
ics can also generate so me turbulence thro ugh spiral way e 
shocks (|Robertslll969[) . rotational s hear (|Schavd 1 2004). 
swing amplified shear instabilities (jHuber fc Pfennigerl 



l200ll ; IWada. Meurer. fc Normanl 12002ft. gravitationa l 
potential energy of the arms (|Elmegreen et al] 120031) . 
and non- linear development of gravi tational instability 
(|Li. Mac Low, fc Klesser] I2005L 12001 . Magnetized disk 
instabilities includ in g the m agnet or otational instability 
dBalbus fc Hawlevl Tl99lL 119981: fSellwood fc Balbud 



Dziourkcvitc h. Elstner. Ri idige T~ 120041 : 
Piontek fc Ostrikerl l2005h , an d the magneto- Jeans 



19991: 



instability ti Wada fc Kodal 12004 iKim fc Ostriken 120061 ; 
iKim et al] I2006D . Finally, there co uld also be a contri- 
butio n from thermal ins t abilit y jKoyama & In utsukal 
2001 iKritsuk fc Normanl 120021 IDib fc Burked 120051 : 



Hcn nebelle fc Auditl 12007) . 



In this paper we attempt a direct observational identi- 
fication of the energy sources that drive the turbulence. 
Our results suggest that the energy input from young 
stars determines the H I line width within the star form- 
ing region of galactic disks, while at large radii, where 
star formation is absent or f eeble, the line width is e i- 
ther sustained by the MRI (|Sellwood fc Balbusl [l999h . 
or by thermal broadening from UV heating. We exam- 
ine whether SN feedback, gas temperature, and MRI can 
maintain effectively the observed ISM line width, assum- 
ing that the energy loss occurs on th e self-dissipatio n 
timescales predicted for th e turbulence (|Mac Lowlll999f h 
IMac Low fc Klessenl (j2004) argue that SNe, among other 
mechanisms of energy injection, are the most effective in 
actively star forming regions, while in the outermost re- 



gions of galaxy disks, where the SFR is low or zero, the 
MRI becomes important in sustaining both the magnetic 
field, even in low density regions, a nd the gas velocity 
dispersion (|Sellwood fc Balbuslfl999l) . 

We use multi-wavelength data in the radio, infrared 
(IR), and far UV (§EJ), to compare detailed maps of the 
ISM kinetic energy density implied by 21 cm line widths 
to maps of the SFR density for a sample of nearby spi- 
ral and dwarf galaxies. In this comparison, we take the 
SFR density as a proxy for the energy input rate to the 
ISM. Relying on the measurement of the velocity dis- 
persion as a tracer of turbulence in the neutral gas, we 
look at both the correlation between SFR and turbulent 
motions in the HI on a pixel- by-pixel basis and at their 
mean values as function of galactocentric distance (§[3]). 
We find that the azimuthal average of the HI velocity 
dispersion decreases radially much more slowly than the 
SFR (§ [4j; specifically, the radial profile of HI veloc- 
ity dispersion approaches the HI thermal value for the 
warm neutral medium where the radial profile of SFR 
truncates. Therefore, we investigate (1) pixel-by-pixel 
how the kinetic energy of the gas turbulence correlates 
with the SFR, especially in terms of SN rate, and (2) 
what role thermal broadening and MRI play in regions 
of low SFR and considerable velocity dispersion (§ [5|. 
In particular, we compare the current turbulent kinetic 
energy of the gas turbulence to the energy input from 
both SNe and MRI (and the efficiencies of this input), 
on the basis that the turbulence dissipation timescales, 
provided by extensive theoretical modeling, approximate 
the energy loss timescales. 

2. DATA 

The main observational data for our analysis are maps 
of the ISM kinetic energy surface density derived from 
high resolution HI data. For a sub-sample we include 
carbon monoxide (CO) emission maps to explore the con- 
tribution of the molecular gas to the kinetic energy bud- 
get. These maps are then compared, pixel by pixel, to 
analogous maps of the SFR surface density, derived from 
UV and thermal IR maps for the same objects. 

We use maps p roduced by The H I Nearby Galaxy 
Survey (THINGS; iWalter et al. l l2008h . a survey of the 
21 cm emission line carried out with the NRAO 6 Very 
Large Array, to obtain the column density £hi and ve- 
locity dispersion ohi of neutral hydrogen. The THINGS 
data have an angular resolution of ~ 7", and a veloc- 
ity resolution of either 2.6 or 5.2 km s" 1 (see Table [lj, 
and a 3 a sensitivity corresponding to a column density 
of AT (HI) > 3.2 x 10 20 cm- 2 . To avoid projection ef- 
fects in the velocity dispersion maps, we consider only 
thos e galaxies that ar e more face-on than ~ 50° (see, 
e.g.. lLerov et al.ll2008l ). leaving 11 disk galaxies in total, 
including three dwarf galaxies (Holmberg II, IC 2574, 
and NGC 4214), listed in Table Q] 

We complement the analysis of the HI gas with the 
mass surface density and velocity dispersion of molec- 
ular hydrogen (H2) gas for about half of our galax- 
ies: NGC 628, NGC 3351, NGC 3184, NGC 4214, 
NGC 4736, and NGC 6946. The H 2 maps are de- 
rived from CO J = 2 — > 1 emission maps, obtained 

6 The National Radio Astronomy Observatory is a facility of 
the National Science Foundation operated under cooperative agree- 
ment by Associated Universities, Inc. 
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at the IRAM 30-m telescope by iLerov et all (|2009h us- 
ing the HERA focal plane array. These maps have res- 
olution of ~ 11" and 2.6 km s" 1 , and sensitivity of 
Eh2 > 4 M© pc~ 2 . To derive Sh2 we assume a con- 
version factor Xqo — 2 x 10 20 cm -2 K km s _1 and a line 
ratio N{J = 2 -> l)/N(J = 1 -> 0) = 0.8. 

The SFR density maps were derived by ILerov et"aT1 
( 2008) from 24 ^m and far-UV (FUV) maps, combining 
dust-enshrouded star formation activity and unobscured 
star fo rmation, respectiv ely. The MIPS 24 /zm emission 
maps ()Rieke et al.l 120041 ) are taken fro m the Spitzer In- 
frare d Nearby Galaxies Survey (SINGS; iKennicutt et al.1 
HSU), and the FUV emission maps are taken from 
the Galaxy Ev olution Explorer Nearb y Galaxy Survey 
(GALEX-NGS; iGil de Paz et~ail[2007t ), for a set of disk 
galaxies in common with the THINGS sample. The 
MIPS instrument has a resolution of ~ 6" at 24 /j,m, 
and a wide areal coverage, which typically covers ~ 2 r25- 
The GALEX FUV images have a resolution (FWHM) of 
5.6" within the A = 1350 — 1750 A band and have a field 
of view of ~ 1.25°. 

3. DATA ANALYSIS 

As a first step in comparing the kinetic energy of the 
gas to possible drivers of turbulence, we convert the pre- 
viously described maps (§ EJ) of HI and H2 gas column 
density and kinematics, and the 24 /im + FUV surface 
brightness into physical quantities. Specifically, we ob- 
tain mass surface density, Shi, Sh2, velocity dispersion, 
chi, ch2, and SFR surface density, Ssfr- Afterwards, 
we proceed with a comparison between the turbulence ki- 
netic energy, obtained from the gas mass surface density 
and velocity dispersion, and the SFR on a pixel-by-pixel 
basis. 

3.1. Mass Surface Density and Velocity Dispersion 

We determine the HI mass surface density, Shi, and 
the velocity dispersion, ohi, by taking the moments of 
fully reduced data cubes for our 11 sample galaxies (the 
prior reduction steps of are described in iWalter et al.1 
2008). We obtain the pixel-by-pixel maps of T,m(x,y), 
v(x, y), and a(x, y), by calculating integrals over the data 
cubes for all the positions on the sky Si(x, y, Vi). All the 
sums are calculated over those channels with signal iden- 
tified by the THINGS masking ([Walter et al.ll2008fl ; we 
only calculated moments for those (x, y) line-of-sight po- 
sitions with signal in at least 3 channels. 

The zero-th moment, Shi, is calculated by integrating 
the H I data cube along the velocity dimension for a total 
number N of velocity channels using: 



-HI 



N 

£ 

i=l 



Si 



(1) 



where Si denotes the signal within the i-th velocity chan- 
nel. Afterwards, Shi is corrected for inclination; the 
adopted inclinations are listed in Table [TJ For each sam- 
ple galaxy, the spectral dimension of the H I data cubes 
is sufficiently broad to encompass all relevant velocity 
values. 

The first moment map, the intensity-weighted mean 
velocity along the line of sight, is given by: 



where Vi is the velocity value of the i-th channel. We 
only use v in this context to calculate the central 2nd- 
moment. 

The variance, corresponding to the square of the line- 
of-sight velocity dispersion, is given by the intensity- 
weighted mean deviation: 
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Throughout our analysis we restrict ourselves to the 
moment definition of Eq. [3J to define the variance. We 
presume that a (Eq. [3]) is the most sensible tracer of tur- 
bulent or more precisely, disordered motions, since it is 
applicable to any line profile. We address the discussion 
of this assumption in § 13.21 

We complement the H I maps with H2 column density 
and kinematics data obtained from CO emission (§[2]). 
Following the same procedure previously adopted to cal- 
culate the relevant HI maps, we obtain the mass surface 
density, £h2 , and the velocity dispersion, cth2 , for the H2 
gas. As described in the following sections, while H2 is 
an important contributor to the total mass surface den- 
sity, its contribution to the total turbulent kinetic energy 
of gas turns out to be rather small in our sample. 

3.2. Realities of Measuring the Velocity Dispersion in 

THINGS 

In the case of perfect signal, the second moment is 
exactly the kinetic energy of HI gas. Even at limited 
signal-to-noise ratio (S/N) but with normally-distributed 
noise, we expect the measured variance to scatter about 
the true value in a well-behaved way, so that by averaging 
measurements over enough area, we can approximate the 
true value. 

As an alternative approach, assuming and then fitting 
a fixed line profile shape (e.g., a Gaussian) fails in two 
ways. First, towards the centers of galaxies, line profiles 
tend to be broad and often contain multiple emission 
components. Second, in the outskirts of galaxies, where 
individual lines of sight tend to have relatively low S/N, a 
simple Gaussian fit tends to underestimate the line width 
by latching onto only the peak of the emission. In the 
low S/N regime, the line profiles are approximately Gaus- 
sian. For a few sample galaxies (NGC 3184, NGC 5194, 
NGC 6946, and NGC 7793) we find that the Gaussian 
a value is systematically lower than the second moment. 
For these galaxies, however, the difference between Gaus- 
sian a and 2nd-moment is lower than ^-2 — 3 km s _1 in 
average. Since we do not know a priori the H I gas line-of- 
sight velocity distribution, and the profiles are not nec- 
essarily Gaussian, the best tracer of kinetic energy, as a 
model independent estimate, is the variance. 

Therefore, we avoid fitting and attempt to minimize 
systematics by measuring the second moment only on 
regions that have been identified to contain signal. The 
production o f thes e "blanked cubes" is described in 
IWalter et al.1 (|2008l ). and is briefly summarized here. 
First, the cubes were smoothed to 30" angular resolu- 
tion. Then, all areas without S/N > 2 in 2 consecu- 
tive velocity channels are blanked (set to 0). This does 
an excellent job of identifying all areas with significant 
emission but still left some false positives. These were re- 
moved by a by-eye inspection that focused on agreement 
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with the overall velocity structure and morphology of the 
galaxy. We test the systematics induced by such blank- 
ing on a face-on galaxy by measuring a on the blanked 
cube and then expanding the blanking mask and remea- 
suring a. We find that small expansions of the mask 
on both directions along the velocity dimension, e.g. by 
±15 km s _1 , have almost no impact on the measured 
dispersion. Larger expansion begin to show systematic 
effects that we want to avoid (particularly at low S/N). 

These blanked cubes minimize the effect of artifacts 
in the THINGS cubes and yield a more robust estimate 
of kinetic energy than fitting an assumed line profile. 
Therefore, the profiles of a and Ef. represent the best 
possible estimates using the available data. However, we 
emphasize that a rigorous attempt to match other obser- 
vations or simulations to our measurements should bear 
in mind (and ideally duplicate or simulate) this blanking 
procedure. 

3.3. Turbulent Kinetic Energy 

With the moment definition of Shi and <thi, the HI ki- 
netic energy per unit area is given by Ek = 3/2 Shi <7hi> 
where the factor 3/2 takes into account all three veloc- 
ity components assuming that the velocity dispersion is 
isotropic. For the galaxies with molecular gas maps, we 
define the total kinetic energy as Ek = 3/2 x (Shi c hi + 
Sh2 c H2 ). I n § HI we describe the relationship between 
the HI (or HI + H2) gas kinetic energy Ek and the Ssfr 
maps in a pixel- by-pixel scatter plot. 

3.4. Star Formation Rate 

We use maps of SFR surface density, Ssfr, derived 
bv iLerov et al.l (|2008l ). from combining the 24 fmx and 
the FUV emission maps ta ken from the SING S and the 
GALEX-NGS surveys (§[2J). ILerov et all (p008h have cal- 
ibrated the SFR represented by the UV and IR emis- 
sion by comparing these maps to Ha, 24 fim,, and Pact 
emission from HII reg i ons a nd young compact stellar 
clusters (jCalzetti et al.ll2007t) . The 24 /an band emis- 
sion is mostly radiation from hot dust heated by the UV 
light from young massive stars and therefore traces dust- 
enshrouded, ongoing star format i on over a time scale 
3 - 10 Mvr (ICalzetti etaTJ 120051: iPerez-Gonzalez et al.1 
2006; iTamburro et al.ll2008l ). although part of the emis- 
sion proceeds from outside the HII regions, thus trac- 
ing older (> 10 Myr) stellar populations. FUV emis- 
sion is mostly photospheric emission from O and B 
stars. It thus complements the 24 /im emission in regions 
poor in dust content, probing therefore low-metallicity 
and older region s of star formation over timescales o f 
t - 10-100 Myr ijOalzetti et al.ll200alSaTim et al.ll200l . 

The calibration provides SFR estimates with an uncer- 
tainty of 40% at most at high SFR, depending on vari- 
ations in geo metry, dust temp erature, and age of stellar 
populations (|Lerov et al.ll2008t ). At low SFR, a substan- 
tial part of the 24 fim. emission proceeds from diffuse 
dust in the ISM, which may not be associated directly 
to recent star formation, with the dust b eing heated 
by ne arby (young and old) star clusters. ILerov et al.l 
(2008) argue that below a fiducial threshold of Ssfr = 
10~ 10 M© yr _1 pc~ 2 the SFR maps represent upper lim- 
its to the true SFR because of the contribution of this 
diffuse dust component. 



4. RESULTS 

In this section we describe our two main results: (1) 
the observation that a radial decline of the HI veloc- 
ity dispersion is pervasive throughout the entire sample 
with little dependence on galaxy type; and (2) that we 
find correlation between the kinetic energy of gas and 
SFR with a slope close to unity and a similar propor- 
tionality constant in all objects. Both results depend on 
the exceptional quality of the data in terms of spatial 
and velocity resolution and the wide field of view. We 
reserve the broader interpretation of our results for § [5l 

In the following, we use units of ^25, whi c h pro - 
vides a conveni e nt no rmalization. ILerov et al.l (|2008h : 
ITamburro et al . (2008) measured the exponential scale 
lengths of near-IR emission from our sample, a good 
proxy for stellar mass, is typically 4.6 ± 0.8 times 
the near-IR scale length, and the SFR maps yield com- 
parable scale lengths to the near-IR. When we normalize 
by ^25 then, it is roughly equivalent to normalizing by 
the scale length of the disk. 

4.1. Radial Profiles of am 

To start, we determine the azimuthally averaged ra- 
dial profiles of HI velocity dispersion, SFR, and gas ki- 
netic energy, by calculating the average values of <thi 
and SFR within annuli of ~ 15" width. The resulting 
chi(»") and Ssfr(^) profiles are shown in Fig.[TJ They ex- 
hibit a radial decline of H I velocity dispersion as a com- 
mon characteristic for all sample galaxies independent 
of their dynamical mass. While previous studies have 
reported this individually for a few disk galaxies (i.e. , 
iBoulanger fc Viallefond|[l992t iPetric fc Rupenll2007D . we 
show for the first time that this is true for a significant 
sample of dwarf and normal spiral galaxies. 

To characterize the ohi radial gradients, we fit a lin- 
ear relation to the radial profiles of <jhi for all r > r25, 
a regime where <j(r) is approximately linear with r for 
all the galaxies of the sample. With a(r) going from 
> 20 km s _1 down to ~ 5 ± 2 km s _1 near the outermost 
observed radius (Fig.Q}, the velocity dispersion decreases 
with radius by ~ 3-5 km s _1 per Ar2s; for comparison, 
azimuthal variations at a fixed radius are ~ 5 km s _1 . 
In Fig. [2j we display and summarize the intercept value 
of <thi at r25 and the H I mass weighted median of cthi 
for all sample galaxies individually. Our observations re- 
solve the outward radial cthi decline well, since a typical 
r25 for our sample galaxies corresponds to ~ 30-40 times 
the 7" HI resolution limit and the typical radial extent 
of the H I emission is 2 — 4 x r25 . 

Remarkably, we find that all galaxies have the same 
HI velocity dispersion at their respective ^5: a general 
value of cr(r25) ~ 10 ± 2 km s _1 , which displays no ap- 
parent trend with the dynamical mass and morpholog- 
ical type, and is consistent with the HI mass- weighted 
median value (a) (see Fig. [5] and Tabled]). The radial 
ohi gradient has in all cases the same sign as the much 
steeper decline of the mean SFR as function of radius as 
shown for comparison in logarithmic scale in Fig. [TJ 

Since the sample galaxies are more face-on than 50°, 
the systematic increase of cthi towards the center cannot 
be due to increasing beam smearing effects. We analyze 
the combined effect of beam smearing, inclination and 
rotation curve on the velocity dispersion by constructing 
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a sample data cube containing signal in only one veloc- 
ity channel per each (x, y) spatial position (i.e. intrinsic 
cm = 0) and characterized by the same inclination and 
rotation curve of NGC 5055 - the most inclined disk of 
the sample with a fast rotation speed. After convolv- 
ing this sample data cube with a kernel of 7" FWHM 
(the resolution of our HI maps), we calculate the ve- 
locity dispersion using second moments (Eq. [3J as done 
throughout our analysis (§ 13. The resulting line broad- 
ening from beam smearing for a disk with the rotation 
curve and orientation of NGC 5055 is only < 5 km s _1 , 
lower than the observed velocity dispersion in NGC 5055 
by ~ 10-20 km s _1 . The fractional contribution of the 
velocity dispersion from beam smearing to the total ob- 
served velocity dispersion is only 20% in the central re- 
gion (r < 1/4 r2s) and at most 10% at larger galactocen- 
tric radii. 

4.2. The Hi Kinetic Energy Density as a Function of 

Radius 

The HI kinetic energy density, Ek, in each pixel ex- 
hibits a clear radial decline, as shown in the full pixel- 
by-pixel distribution (Fig. [SJ, where the black contours 
and color scale indicate the density of pixels at each H I 
kinetic energy and radius, while the red contours show 
the sum of atomic and molecular kinetic energy. The 
latter diverges from the HI kinetic energy in the inner 
parts of some galaxies, indicating that the cold molecular 
gas contributes to the kinetic energy budget. Since every 
galaxy of the sample includes ~2x 10 5 pixels with sig- 
nificant signal, we display in Fig. [3] the density contours 
of the data points. 

We include the analysis of the H2 mass surface density 
and velocity dispersion derived from the CO emission 
for a few sample galaxies (§ 0), to quantify the contri- 
bution of the molecular gas to the total kinetic energy 
at high H2-to-HI mass ratio. The total kinetic energy 
Ek = Em + En2 is plotted in Fig. [3] (with red con- 
tours) for the galaxies with CO data. For the galaxies 
NGC 4736, NGC 5055, and NGC 6946, E m is compa- 
rable to Em or even dominant in the central regions of 
galaxy disks; for these galaxies the H2 mass and spa- 
tial extent is considerable. For those galaxies where the 
HI gas dominates the total gas mass, i.e., NGC 628, 
NGC 3184, NGC 3351, and the dwarf galaxies Holm- 
berg II, IC 2474, and NGC 4214, characterized by little 
or no detected molecular gas, the molecular gas does not 
contribute much to the total kinetic energy. 

4.3. Correlation between Ek and Esfr 

In Fig. [5] we compare pixel-by-pixel the relation be- 
tween the kinetic energy density of the HI gas, Ek = 
3/2£hiOhi> an d the SFR surface density - a proxy for 
the energy input rate by SN. We find that in all galax- 
ies these quantities are well correlated with a slope close 
to unity and with no evident dependence on dynamical 
mass (Fig. [4]). Note that a considerable fraction of all 
data points in Fig. 0] lie below the noise estimated for 
the Ssfr maps (§ 13 ,4|) . These data points lie in the out- 
ermost parts of galaxy disks (r > 2 x ^5), where there 
is little or no ongoing star formation (cf. Fig. [3]). 

We note that the slope of all correlations in Fig. 2] flat- 
tens at high Ek and £sfr, although we argue in § l4.2l that 
this is not caused by a higher abundance of molecular gas 



at higher Ssfr- The total kinetic energy Ek — £'hi+-E'h2 
is plotted in Fig. 2] (with red contours) for the six sam- 
ple galaxies with H2 data. Only for NGC 5055 and 
NGC 6946 is the kinetic energy of the molecular gas, 
Eu2 important; in those cases the Em + £h2 vs. Esfr 
relation is linear and the slope is close to unity. 

5. DISCUSSION 

What scenario does our data support for producing the 
observed line widths? SN-driven turbulence seems likely 
to be the dominant factor broadening line widths within 
the radius of active star formation, since we find that the 
level of predicted SN energy is sufficient to account for 
the turbulent kinetic energy implied by the line width as 
a function of radius. 

The radial slopes of SN energy and kinetic energy of 
the neutral (HI and H 2 ) gas agree qualitatively so that 
the kinetic energy of the gas is proportional to the lo- 
cal star formation rate. Yet, the fact that the HI ve- 
locity dispersion approaches its thermal value of roughly 
6 km s _1 well beyond the radius of detectable star forma- 
tion indicates that either (1) the line broadening is due to 
UV heating, with a warm neutral medium temperature 
of ~ 5000 K resulting in cthi ~ 6 km s _1 , or (2) the gas 
is actually turbulent and another mechanism such as the 
MRI is driving the turbulence. We now explore whether 
our new data support this scenario. 

In the following, we compare the observed line widths 
with the most plausible mechanisms for generating them: 
SNe fSO]). UV heating (§ E2J), and MRI (§p. SNe 
and MRI both produce broad line widths by driving tur- 
bulence, while UV heating can produce thermal broad- 
ening. The required energy injection rate depends both 
on the kinetic or thermal energy of the gas, which can be 
derived from the observed line widths and the turbulence 
decay or cooling timescales, which must be derived from 
models. More precisely, the gas kinetic energy implied 
by the linewidth consists of a combination of turbulence 
and the thermal energy associated with the (warm) gas 
temperature, i.e. Ek = E tur b + ^thcrm- If the thermal 
broadening is much less effective than the turbulence, 
then E k ~ -Eturb- 

If the gas turbulence is mainly driven by SNe and MRI, 
then we expect the dissipation rate of turbulence to equal 
the sum of the energy input rates of SNe and MRI, 

Ek — £sn -Esn/tsn + emri^mri/tmri, (4) 

where £sn and £mri are the efficiencies, and tsn and 
tmri are the decay times of turbulence driven by the two 
mechanisms. Different mechanisms can result in different 
decay rates b ecause they have different driving scales and 
magnitudes (|Stone et al.lll998t iMac Lowj|1999l ). 

5.1. Supernova Energy 

Assuming steady state equilibrium between the energy 
input rate from SNe to turbulent gas motions and the 
energy loss rate from dissipation of this turbulence, then 
the resulting kinetic energy Ek — esN-EsNTSJV , where 
Esn is the rate of released SN energy, which we estimate 
from the SFR, and the SN feedback efficiency 6sn is the 
fraction of S N energy c o nvert ed to turbulent motions in 
the cold gas. IMac Low! (|1999f ) finds that the dissipation 
rate of turbulence depends on the driving scale A and the 
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velocity dispersion a as 



td — 9.8 (Aioo/cio) Myr, 



(5) 



where Aioo = A/100 pc and aio = cr/10 km s . Numer- 
ical si mulations of SN-driyen turbulence yield A = 100 ± 
30 pc (| Joung fe Mac Low|[200l lAvillez fc BreitschwerdtJ 
l2007h . and our own analysis gives an average velocity 
dispersion er = 10 km s _1 (§ 14. ip . The SN energy input 

rate, -Esn, can be estimated from the SN rate implied by 
our SFR maps. The SN rate per unit area, T], depends 
on the fraction /*^sn of all recently formed stars that 
terminate in core-collapse SNe: 



SFR 

v = -i—c x /* 

(m) 



*SN, 



(6) 



where (m) is the average mass of stars of the popu- 
lation. We assume that only those stars in the mass 
range (8-120) M Q can form core-collapse SNe. The SFR 
maps used in our analysis assume an initial mass func- 
tion (IMF) (f)(m) = m~ a , where a = 1.3 for the mass 
range (0.1 — 0.5) M q and a = 2.3 for the mass range 
(0.5 - 120) M (cf. iLeitherer et~aTlll999t ICalzetti et all 
2007). Then, the SN fraction 



JsTq" 8 <t>( m ) dm JoTm® mcf>(m) dm 

JoTmI «K m ) dlJl JoSM^ <H TO ) dm 



-120 M f 



+SN 



/(m) 



(7) 

yielding f*^sN/(m) — 1-3 x 10 2 M Q . If we were to re- 
duce the upper mass limit to 50 M Q yields /*_ + sn/(™) — 
1.2 x 10 -2 Mq 1 , and an upper mass limit of 20 M© yields 
0.9 x 10 -2 Mq 1 . The effect of these variations on r\ is 
unclear, because the upper mass limit of the IMF also 
affects our translation of UV and IR light into SFR. The 
UV and IR maps are primarily sensitive to high mass 
stars. Therefore for a lower upper mass limit, therefore, 
less UV and IR light is emitted per unit star formed and 
if the upper mass limit is actually lower than we have 
assumed, then we have underestimated the true SFR. In 
calculating 77, these higher SFR and lower /„^sn/(to) 
have opposite effects, leaving the impact of changing the 
upper mass limit on 77 unclear. As the upper mass limit 
decreases, we estimate less intrinsic UV and IR emission 
per high-mass star, which we have used to construct the 
maps. We neglect the contribution of type la SNe whose 
rate is ~ 1/3 of the core-collapse rate for the morpholog - 
ical types of our sample galaxies (|Mannucci et al.l f2005). 
In these circumstances, the SN rate can be straightfor- 
wardly calculated as a function of SFR, 77 = 7/(Esfr), 
and, depending on the adopted assumptions, the SN rate 
is characterized by an uncertainty of a factor ~ 1/3. We 
assume that for each SN explosion only a fraction esN < 1 
of 10 5 1 erg - rough ly the energy released by a single SN 
event (Heiles 1987) - is converted into turbulence. Then 
in steady state the kinetic energy of the gas turbulence 
E k =rjx (e SN 10 51 erg)r D . 

In Fig.s [3] and |4j we compare our estimate for the to- 
tal energy from SNe produced within a single turbulent 
decay time td = 9.8 Myr with Ek, the kinetic energy 
derived from the line width, as a function of radius and 
of SFR, respectively, for different values of SN efficiency. 
Fig. [3] shows the azimuthally averaged SN energy de- 
creasing as a function of radius. A universal £sn would 



produce a linear correlation between Ek and Esfr as 
shown in Fig. 01 where lines of unity slope and constant 
values of 6gN represent SN energy input as a function of 
SFR. Fig. [3] and Fig. [4] show that, at least in the high 
star formation regime (Ssfr > 10~ 9 M Q yr _1 pc -2 ), 
the data points typically agree with an efficiency 0.1 < 
e SN x (10 7 yr/ TD ) < 1. 

This estimate of the efficiency is consistent with 
numerical simulation s that estimate (csn) — 0.1 
(|Thornton et al.lll998fl . In the low SFR regime (£ S fr < 
10~ 9 M Q yr _1 pc -2 ), many data points in Fig.s [3] and [4] 
lie close to the line of maximum efficiency esn = 1- At 
^SFR < 10 -10 M yr _1 pc~ 2 , in particular at large 
galactocentric radii, the estimated Ssfr rates fall to 
within the noise (see § I3.4p . whilst the measured Ek is 
typically well constrained. Nevertheless, esn = 1 remains 
problematic, since it would imply that all the kinetic en- 
ergy of SN remnants would be deposited as kine tic en- 
ergy of the gas. In fact iTenorio-Tagle et al.l (|1991h argue 
that the expected SN efficiency esN should be at most 
~ 0.5. Values graeter than 0.5, therefore, either imply 
that other sources inject energy into the ISM (see fol- 
lowing sections), or that the dissipation timescales are 
shorter than we assume. In the next sections, we exam- 
ine alternative mechanisms to explain the line width in 
these regions. 

5.2. Thermal Broadening 

The temperature of the warm phase of the HI is 
maintained by UV heating. However, neutral gas never 
reaches temperatures high enough to explain line widths 
as high as 10 km s , which ar e instead attr i buted 
to supersonic tur bulent motions (|Wolfire et al.l 120031 : 
iLacour et all 120051 ). As a m > 10 km s -1 for r < r 25 
for all our sample galaxies, we argue that thermal ef- 
fects can be neglected within the star forming radius. In 
regions of active star formation, stellar winds and ion- 
izing radiation from stars appear less effective together 
than the S Ne from the same stellar po pulation at driving 
turbulence (|Mac Low fc Klessen|[200l . Thus, SN driven 
turbulence looks likely to dominate there. 

In the outer regions of H I disks, on the other hand, the 
velocity dispersion approaches its thermal value. There, 
the warm neutral atomic phase temperature is typically 
~ 5000 K, as measured , e.g., in the solar neig hborhood 
(|Heiles fc Trolandll2003l : iRedfield fc Linsky||2004h . where 
the SFR is a few xlO -9 M Q yr _1 pc -2 . This tempera- 
ture gives a thermal width of cthi ~ 6 km s _1 . Typical 
temperatures of the warm medium can eve n be as high as 
~ 85 00 K (corresponding to ~ 8 km s _1 ; IWolfire et al.l 
1995). If thermal broadening is effective, the observed 
cm ~ 6 km s _1 outside the star forming radius does not 
necessarily involve turbulence. However, such tempera- 
ture levels, especially at r > 2 x r 2 5, where the SFR is 
below 10~ 10 M yr _1 pc -2 (cf. Fig. QJ, still require a 
continuous UV background source warming up the H I 
disks of galaxies. At r > 7-25, the source of such UV radi- 
ation is presumably not local; it could be extragalactic. 

The local thermal pressure estimate and the actual 
velocity dispersion in the outer parts of HI disks are 
consistent with the existence of a wa rm phase for the 
gas. For example, (jLeroy et al.l 120081 ) estimate a local 
pressure P/k ~ 300 K cm -3 at r ~ 2 X r 2 5 for the 
galaxy NGC 4214, assuming hydrostatic equilibrium (cf. 
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Elmegrecn 1989). At a radius 2 x r25, the HI velocity 
dispersion in NGC 4214 is ~ 7 km s" 1 , corresponding to 
a temperature of ~ 5900 K, which, solving P oc no 2 , 
yields a density n ~ 0.05 cm -3 . This is consistent 
with the gas i n a w arm phase according to Figure 7 in 
IWolfire et all (|2003h . where a pressure of ~ 300 K cm 
corresponds to n ~ 0.04 cm -3 , at least for the outer parts 
(r ~ 15 — 18 kpc) of the Milky Way, which likely have a 
similar UV background. For pressure values lower than 
300 K cm -3 , all the gas is warm. 

5.3. MRI Energy 

If external UV heating is insufficient to maintain a 
warm phase, then a non-stellar energy source is needed 
to explain the observed turbulence (§ I4.2|) . MRI is a 
plausible candidate. It develops in differentially rotat- 
ing disks with angular velocity decreasing outwards, as 
in all but the smallest galactic disks, as long as some 
weak magnetic field threads the disk. It has been argued 
to sustain both the interstellar gas turbulence and the 
galactic ma gnetic field to a few /xG (|Piontek fc Ostrikerl 
120051 120071 ) . MRI requires a minimum magnetic field as 
low as 10~ 25 G to originate, w hich is much lower than 
the s eed galactic magnetic fields dKitchatinov fc Riidigerl 
12001 . 

Following the same line of reasoning as § 15. 1[ if 
we assume that a fraction emri < 1 of the MRI 
energy is transformed into kinetic energy, then in 
steady state the observed kinetic energy must equal 
the MRI energy input within the turbulence decay 
time: Ek = emr.i Embi Tmri- Theoretical calcu- 
lations dSellwood fc Balbud 119991 : iMac Low fc Klessenl 
120041 : IMac Low! l2008| ) estimate the production energy 
rate of MRI to be 
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(8) 



where 51 = v/r is the angular velocity, h z is the vertical 
thickness of the HI disk, and B is the magnetic field. In 
our analysis, we assume a constant thickness h z = 100 pc 
for all galaxies. This is ap propriate for the in ner Milky 
Way out to the solar circle ijWolfire et al .1 120031 , although 
in the outer regions of disks a nd in dwarf galaxies h z 
may be higher (up to ~ 300 pc: IWalter fc Brinkdll999l ). 
We also assume a constant magnetic field B = 6 /iG 
as a typical galactic m agnetic field (|Beck et al.l 119961 : 
iHeiles fc Trolandl I2005T ). Taking the turbulent decay 
timescale again to be tmri — 9.8 (Aioo/fio) Myr we can 
estimate the driving scale A = A c . The critical wave- 
lengt h for the fastest MRI growth (|Balbus fc Hawlevl 
1998) 



A c = 2tt va 
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(9) 



where a = d\nv/d\nr, and the Alfven velocity v\ = 
B 2 /4wp. For a typical galactic magnetic field of B — 
6 [iG and a density of 2 x 10~ 24 g cm~ 3 , A c ~ 10 2 pc. 

In Fig. [3j we compare the energy produced by MRI 
within a decay time of turbulence, £mri -E-mri tmri, with 
the observed kinetic energy Ek as a function of radius 
for different values of the MRI efficiency. Note that al- 
though both £mri and tmri are functions of radius, the 



MRI energy input only decreases with radius slowly, not 
exponentially as Ek- In Fig. [3l we show a comparison 
between the SN energy f£j |5.1[) and the MRI energy plot- 
ted as a function of radius, and indicated with green and 
blue solid lines, respectively. We also plot in the pixel- 
by-pixel Ek vs. Ssfr plot of Figure |4] the average value 
of the MRI for pixels with each Esfr, for £mri = 1. 

While SNe can be identified as the dominant source of 
energy at Esfr > ~9 M Q yr _1 pc~ 2 , the MRI contri- 
bution becomes important at Ssfr < —9 M a yr _1 pc~ 2 
and dominant at large galactocentric radii. On the other 
hand, the kinetic energy of turbulence within r25 is much 
higher than the predicted MRI energy, indicating that 
MRI can not account for the observed gas turbulence in 
regions of active star formation. 

The observed Ek is higher than the SN energy for 
the dwarf galaxies Holmberg II and, more significantly, 
IC 2574. In these two c ases, the gas turbulence cannot 
be explained by SNe (cf. ISt ammirovic & Lazarian1l2001t 
iDib fc Burkertjl2005HPasauali et al.ll200l : still, the MRI 
could account for the observed regime of turbulence. Our 
analysis suggests that MRI could dominate regions of low 
SFR, and rules out MRI as an effective turbulence driv- 
ing mechanism at high SFR. 

5.4. Robustness of the Approach 

Three major sources of uncertainty enter our analysis. 
The first and largest of these are the empirical conver- 
sions from UV and IR emission to Ssfr, which may in- 
troduce up to a 40% uncertainty (§ 13. 4p . The conversion 
from Ssfr to SN rate relying on a universal IMF intro- 
duces an additional 30% uncertainty (§ I5.ip . A second 
uncertainty enters from the reliance on numerical simu- 
lation to evaluate the driving scale of SN-driven turbu- 
lence and thus the dissipation timescale tsn- The simu- 
lation s give a 30% error for the i r estimate of the driving 
scale ([Joung fc Mac Low! [20061 : lAvillez fc Breitschwerdl 
2007). The dissipation timescale also depends on the 
mean velocity dispersion, which is, of course, not con- 
stant as assumed in our estimate. A third uncertainty 
enters in our estimates of MRI energy, where we assume 
a constant magnetic field and vertical thickness of the 
disk for all the galaxies of the sample, which may not 
be the case. The vertical thickness of the gaseous com- 
ponent in galaxies increases outwards and may double 
the value assumed in our analysis of h z — 100 pc. The 
typical strength of magnetic fields observed in the Milky 
Way and other galaxies declines slowly as a function of 
galactocentric distance, although the variations of the 
magne tic field can be much larger azimuthally than ra - 
dially ([Fletcher et al.ll2004t lHan et al.l[200llBecldl2007] ). 
These three sources of errors could in principle explain 
those data points lying near £sn ^ 1 and cmri ^ 1 ■ 

Aside from the potential sources of uncertainties, if we 
were to interpret our empirical finding by taking into ac- 
count only SNe explosions and MRI, still a minor part of 
the observed turbulence would require csn and cmki effi- 
ciencies uncomfortably high, as high as 100%. Therefore, 
we do not exclude that other mechanisms could be effi- 
ciently driving some turbulence in the gas. Potentially, 
within the star forming radius (< ^5) stellar winds could 
be the most effective, while in the outermost regions of 
HI disks, i.e. r > 2 x r-25, where the SNe and the star 
formation effects are not likely to produce feedback, a 
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floor level of HI velocity dispersion of ~ 6 km s _1 could 
be attributed to thermal broadening. Observation of the 
H I in absorption at large galactocentric radii might ulti- 
mately help in determi ning whether the gas i s cold and 
turbulent or warm fsee lDickev fe Brinkdll993h . 

5.5. Is Shi Controlling the Ek vs Ssfr Relation? 

The gas kinetic energy, Ek, is the correct physical 
quantity to study the energy balance in the ISM. How- 
ever, gas surface densit y and SFR are well known to cor- 
relate (e.g. !KennicuttlTi998t iBigiel et al.l 120081 ). so one 
may wonder whether this drives our observed correlation 
between Ek and Ssfr- In other words, the observed co- 
variation between Ek and Ssfr might result from the 
two facts: at higher gas mass density, galaxies form stars 
at higher rate, and higher gas mass bears higher kinetic 
energy. In order to verify that the correlation between 
E/~ and Ssfr ^ s n °t controlled by Shi, we remove the 
effect of Shi on the Ek vs Ssfr relation by calculating 
the partial correlation coefficient: 



P12.3 
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where P12.3 is the partial correlation between x\ = 
log Ssfr and X2 = log-Efc while controlling for 23 = 
log Shi, and pij is the Pearson's correlation coefficient 
between two data sets Xi and Xj. Here, considering the 
quantities in logarithmic scale is convenient as they are 
correlated as power laws. Retaining only data points 
above a fiducia l value for the H I mass density, i.e. 
S H i > 3 M© (cf. IBigiel et~aTll2008| ). we obtain the values 
listed in Table [2] for our sample galaxies. If the corre- 
lation between Ek and Ssfr were completely controlled 
by Shi we would expect P12.3 — 0. Yet, we find that 
0-2 < P12.3 < 0.6 for our sample galaxies, indicating that 
the correlation between Ek and Ssfr is real. Equiva- 
lently, we show that the correlation Ek vs Ssfr at con- 
stant Shi holds a positive slope in Fig[5j in which we 
remove the contribution from Shi to Ek by subtracting 
the average ehi within bins of Shi- The residuals 
Ek — (-Efc)sHi vs Ssfr exhibit a positive correlation, in- 
dicating that it is not Shi alone that determines the ob- 
served Ek vs Ssfr correlation. Positive slopes in Fig. [5] 
and positive partial correlation coefficients imply that 
there is a real, physical relationship between Ek and SFR 
even at fixed Shi- The relatively weak slopes in Fig. [5] 
indicate that higher gas mass density correlates indeed 
with higher kinetic energy, simply because it generates 
more star formation. Fig. [5] and the partial correlation 
coefficients allow us to detect a relationship between Ek 
and SFR that is independent of Shi- Comparing Fig. [4] 
and Fig. [5j however, it is clear that most of the cor- 
relation between Ek and SFR in Fig. [5] closely involves 
Shi (the distributions in Fig. [5] are very flat compared 
to those in Fig. |4j. The basic effect seems to be that 
higher Shi results in higher SFR, which creates more Ek- 
Also, note that if the turbulence, as traced by Ek, were 
effectively suppressing star formation, we would have ob- 
served a negative correlation here. 

5.6. Does Turbulence Drive Stochastic Star Formation? 

The data analysis suggests that the SFR drives the 
HI turbulence through SN feedback. However, the ob- 



served correlation could also be interpreted in the oppo- 
site logical direction, i.e., that the turbulence is driv- 
ing star formation. In f act, as it has been argued 
(|Mac Low fc Klessenir2004l) . turbulence in the ISM has 
a dual role: (1) to quench star formation by providing 
pressure support to the ISM and preventing collapse, and 
(2) to promote star formation by generating stochastic 
super-critical density enhancements. If the turbulence 
were to drive substantial stochastic star formation, we 
would indeed expect a positive correlation between <thi 
and Ssfr- However, Fig. Q] shows that cthi and Ssfr 
occupy quite different dynamic ranges. While the Ssfr 
ranges over several orders of magnitude, <thi ranges from 
~ 20 to ~ 5 km s" 1 and is characterized by large az- 
imuthal variations. Although Ssfr positively correlates 
with <thi on galactic scales, the large azimuthal varia- 
tions imply that Ssfr is not well defined for any given 
value of cthi- This does not preclude turbulent induc- 
tion of star formation in individual regions, but does 
suggest that this process does not dominate over large 
scales. The physical explanation might be that supercrit- 
ical density fluctuations are often dispersed on timescales 
shorter than the free-fall time, arresting the collapse 
dKlessen et alj|2000t lElmegreeril200l I Joung fc Mac Low! 
12006ft . 

6. EFFECTS OF SPIRAL ARM KINEMATICS AND TIDAL 
INTERACTIONS 

In the following we discuss other possible mechanisms 
to produce ISM turbulence, such as spiral arm kinemat- 
ics, tidal interactions, and streamers. 

Table [TJ which lists the morphological types of the 
galaxies of our sample, shows that there is no evident 
trend of the mean velocity dispersion, (a), among in- 
dividual galaxies or morphological type. Spiral galaxies 
with strong spiral pattern, e.g. NGC 628 and NGC 3184, 
have similar values of typical HI velocity dispersion as 
galaxies with no clear spiral structure, e.g. Holmberg II 
and IC 2574. Since the spiral arm strength should vary 
within the sample, we argue that the spiral arm kinemat- 
ics in our sample galaxies are not an important effect in 
driving turbulence into the ISM. 

Although the galaxies Holmberg II and IC 2574 belong 
to the M81 group, they do not show signatures of tidal 
distortion. In our sample, only NGC 5194 is an interact- 
ing galaxy. The HI velocity dispersion in NGC 5194 is 
significantly higher than the average for the galaxies in 
the sample (see Fig. [2]). On the basis of our results, we 
speculate that the tidal interaction with the companion 
NGC 5195 enhanced the SFR in the disk of NGC 5194, 
which has consequently driven the velocity dispersion in 
the HI gas to higher values. 

Extended streamers characterize the galaxies 
NCG 4736 and NGC 5055. Their radial profiles of 
the HI velocity dispersion exhibit a local increase 
outside the radius of active star formation, i.e., at r ~ 4' 
(~ r 25 ) for NCG 4736 and at r ~ 9' (~ 1.5 r 25 ) for 
NGC 5055. However, these local peaks in cthi(^) and 
the streamers have different galactocentric locations, 
corresponding to r ~ 8' and r > 11' for NCG 4736 
and NGC 5055, respectively. Therefore, we argue that 
the presence of extended streamers is not likely to be 
connected to higher HI velocity dispersion. 

7. CONCLUSIONS 
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Combining high quality maps of HI column density 
and line width provided by THINGS for a sample of 
dwarf and spiral galaxies, we obtain the following results. 

1. The HI velocity dispersion, <thi, declines uniformly 
as a function of galactocentric distance in all ana- 
lyzed galaxies. 

2. At r25, the edge of the star- forming region, the HI 
velocity dispersion u(r25) — 10 ± 2 km s _1 , which 
is consistent with the mass- weighted median HI 
velocity dispersion (a). These findings are inde- 
pendent of the dynamical mass of the galaxy and 
of their morphological type. 

3. Within the radius of active star formation (r < 
r25), the estimated SN rate and the correspond- 
ing energy input rate are sufficient to account for 
the bulk of observed kinetic energy of turbulence. 
For those galaxies of the sample with considerable 
H2 gas, the SNe can well account for the com- 
bined HI and H2 gas turbulence. In this region, 
the observed instantaneous kinetic energy of the 
HI gas is consistent with the balance between the 
energy input from the total number of SNe calcu- 
lated from the observed SFR and the turbulent dis- 
sipation predicted by numerical models. The pro- 
portionality between gas Ek and SN energy input 
rate derived from the SFR provides direct evidence 
that HI turbulence comes from SNe in regions of 
active SFR. The resulting SN feedback efficiencies 
are typically £sn x (10 7 yr/775) ~ 0.1 at SFR levels 
Ssfr > 10~ 9 M yr _1 pc~ 2 , with the dissipation 
timcscale of turbulence td — 10 yr. 

4. Within the star forming disk (r < ^5), neither 
thermal broadening nor MRI can produce the ob- 
served HI velocity dispersion. At low SFR, Xsfr < 
10~ 9 M Q yr -1 pc~ 2 , corresponding to large radial 
distances (r > r 2 s), an additional mechanism driv- 



ing the HI velocity dispersion is required to avoid 
SN efficiencies 6sn > 1. 

5. The thermal broadening of the spectral lines, asso- 
ciated to a temperature of ~ 5000 K, may be able 
to explain the observed cthi ~ 6 km s~ 1 in the out- 
ermost regions of H I disks in our sample galaxies, 
if the required UV radiation to maintain these tem- 
peratures is present. The energy input from MRI 
can account for the kinetic energy observed in re- 
gions of low SFR, Ssfr < 10~ 9 Mq yr" 1 pc~ 2 , at 
large galactocentric distances. 

6. We can not unambiguously separate the tem- 
perature of the warm and non-turbulent neutral 
medium from the effect of MRI stirring the ISM. 
Both mechanisms are equivalently plausible drivers 
of the H I velocity dispersion observed in the outer 
parts (r > ^5) of galaxy disks. We suggest that 
testing the HI line profiles of the gas against a 
bright background source could ultimately clarify 
whether the gas in regions of weak star formation 
is uniformly warm, or contains a cold, turbulent 
phase, presumably stirred by MRI. If the gas is ac- 
tually turbulent, the gas kinetic energy for both 
high and low star forming regions is consistent in 
all cases with realistic values of esN and €mri effi- 
ciencies, suggesting that the feedback provided by 
both SN explosions and MRI is sufficient to drive 
the bulk of the observed HI turbulence. 
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Fig. 1. — Radial profiles of Hi velocity dispersion, <jhi, for the galaxies of our sample (§ 14. li . The filled circles represent azimuthal 
averages of cthi calculated from the second moment (Eq. [3J as a function of radius r in units of T25- The error bar on the top right of 
each panel shows the average standard deviation of azimuthal scatter of ohi- The blue crosses denote the radial profiles of log Ssfr with 
the scale i ndic ated on the right side of the plot. Here, only those values above the noise threshold (log Ssprj/pVlp) yr —1 pc~ 2 ] > —10) are 
plotted ( § 13.41 1. The labeled galaxies are sorted by increasing dynamical mass (u ma x dc Blok et al. 2008; Tamburro ct al. 2008). 
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Fig. 2. — Mean values of Hi velocity dispersion for our sample galaxies (§ 14,11 1, Top panel: the fille d cir cles and the corresponding 
error bars represent the value of the Hi velocity dispersion at r25, cr(r25), from the linear fit in Fig HI f§ 14,11 1: the squares represent the 
Spji-weighted median (cr) with the error bars indicating the 68% confidence interval (cf. Table [l}. The solid horizontal line represents the 
median value (10 ± 2 km s _1 ) of all o"(r2s) for the sample galaxies. The labeled galaxies are sorted by increasing dynamical mass (^ ma x), 
Bottom panel: the squares represent the Em-weighted median (a) as in the top panel; the galaxies are sorted by Hubble type (cf. Table [lj. 
The solid horizontal line represents the median of all (a) ~ 10 ± 2 km s — 1 . 
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Fig. 3. — Radial distributions of Hi kinetic energy per unit area, E^ = 3/2 Ehi o"ht 2 , for the galaxies of our sample (§ 14.21 1. The number 
of pixels falling at each position in the plot is represented in color scale and contoured in black. The red contours represent the sum of 
kinetic energies for Hi and H2 gas, which diverge from the Hi contours only near galactic centers. The solid green lines represent the total SN 
energy released over the self-dissipation timescale of turbulence for values of the SN efficiency egN = 1 (above) and 0.1 (below), as discussed 
in § 15. ll The SN energy estimates plotted here retain only those values above the Esfr noise threshold, EsFR. > —10 Mq yr _1 pc~ 2 
( § 13. H . The blue solid lines indicate MRI energy produced over the c orresponding timescale for turbulence dissipation for values of the 
MRI efficiency emri = 1 (above) and 0.1 (below) as discussed in § 15.31 The gray horizontal shadow indicates the noise level of Ey.. 
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Fig. 4. — Pixel-by-pixel scatter plot of = 3/2 Ehi o"hi 2 vs Sgp^ for the galaxies of our sample (§ 14.31 1, The number of pixels at each 
position is represented in color scale and contoured in black. The red contours represent the sum of kinetic energies for Hi and H 2 gas . The 
solid lines of unity slope represent the SN energy input for different SN efficiency valu es es n = 1, 0.1, 0.01 (top to bottom; see §EU- The 
blue solid line indicates the MRI energy input at maximum efficiency emri = 1 (see § I5.3|l . The error bar of a typical pixel is indicated i n 
the upper right panel. The gray vertical and horizontal shadows indicate the noise level of Esfr an d E^, respectively (Lcroy et al. 2008). 
The last panel represents the superposition of data points for all the galaxies of the sample. 
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FlG. 5. — Pixel-by-pixel scatter plot of residuals E^ — (-Efc)EHl vs EgpR after subtracting the average of at constant Shi- The solid 
lines represent the linear fit to data. This procedure is equivalent to looking at the correlation between E^ and EgpR at constant Shi- 
T he re siduals are positively correlated with Egp^, indicating that Ejji does not determine the observed Ej. vs Esfr correlation alone 
(8 15.51 . Beside this graphical illustration, we perform the analytical test of the actual correlation between E^ and Esfr while controlling 
for Shi through calculating the partial correlation coefficent, as described in the text (see also Table l2l for a list of the partial correlation 
coefficients). 
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Obj. name 


b 


t 


»"25 


dv 




0"25 


da/dr 


morph. 


morph. 




(Mpc) 


(°) 


(') 


(km s" 1 ) 


(km s" 1 ) 


(km s _1 ) 


(km s' 1 r" 1 ) 


code 


type 




(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


Dwarf galaxies 


HO 11 


3.39 


41 


3.3 


2.6 


9.4 


9.9 ± 0.5 


-2.5 ± 0.5 


10 


Irr 


IC 2574 


4.02 


53.4 


6.44 


2.6 


9.7 


10.8 ±0.7 


-6.1 ±0.7 


9 


Sm 


NGC 4214 


2.94 


43.7 


3.38 


1.3 


8.9 


9.2 ±0.3 


-2.1 ±0.3 


10 


Irr 








Normal spiral 


galaxies 










NGC 628 


7.3 


7 


4.77 


2.6 


8.0 


7.7 ±0.4 


-1.8 ±0.4 


5 


sv 


NGC 3184 


11.1 


16 


3.62 


2.6 


10.4 


9.7 ±0.9 


-5.2 ±0.9 


6 


Sc 


NGC 3351 


9.33 


41 


3.54 


5.2 


10.6 


9.6 ±0.9 


-2.6 ±0.9 


3 


Sb 


NGC 4736 


4.66 


41.4 


3.88 


5.2 


12.0 


12.1 ± 1.4 


-4.7 ± 1.4 


2 


Sab 


NGC 5055 


7.82 


59 


6.01 


5.2 


12.6 


11.4 ±0.3 


-1.2 ±0.3 


4 


Sbc 


NGC 5194 


7.77 


42 


4.89 


5.2 


17.7 


17.0 ± 1.4 


-3.3 ± 1.4 


4 


Sbc 


NGC 6946 


5.5 


32.6 


5.35 


2.6 


10.1 


8.6 ±0.9 


-1.0 ±0.9 


6 


Sc 


NGC 7793 


3.82 


50 


5.0 


2.6 


11.4 


9.7 ± 1.1 


-2.6 ± 1.1 


7 


Scd 



TABLE 1 ^ 

THINGS TARGET GALAXIES. (1) ADOPTED DISTANCE (IWALTER ET AL.|[2008D ; (2) INCLINATION (IDE BLOK ET AL.||2008I) ; (3) SEMI-MAJOR 
AXIS OF THE 25 MAG ARCSEC -2 ISOPHOTE IN THE B BAND OBTAINED FROM THE LEDA DATABASE (URL: | HTTP : / /LEDA . UNIV-LYON 1 . FR/| ; 
(4) VELOCITY RESOLUTION OF TH E Hi DATA CUBES: (5) Hi MASS- WEIGHTED MEDIAN OF THE Hi VELOCITY DISPERSION; (fa ) MI VELOCITY 
DISPERSION o(r25) AT T25 (§ 14.11 FIGURE CAPTION FlG.0; (7) SLOPES OF RADIAL (JHI PROFILES IN UNITS OF T25 (§ 14.11 1; (8) 
MORPHOLOGICAL CODE FOR THE REVISED DE VAUCOULEURS TYPE (LEDA); (9) MORPHOLOGICAL HUBBLE TYPE (LEDA). 



object 


P12.3 


HO II 


0.54 


NGC 4214 


0.32 


IC 2574 


0.19 


NGC 7793 


0.17 


NGC 4736 


0.58 


NGC 6946 


0.52 


NGC 5055 


0.26 


NGC 3351 


0.62 


NGC 628 


0.52 


NGC 5194 


0.44 


NGC 3184 


0.49 



TABLE 2 

Partial correlation coefficients calculated in order to test the act ual correlation between E k and £g FR while 

CONTROLLING FOR S H I (§ 15.51 ). 



